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Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, 
have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here 
we propose a new way of launching high density and high velocity, plasma jets using multiple intense 
laser beams in a hollow ring formation. We show that such jets provide a more flexible and versa- 
tile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic 
simulations, we demonstrate that the collimated jets can achieve much higher density, temperature 
and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, in- 
stead of focused onto a single spot. We carried out simulations with different ring radii and studied 
their effects on the jet properties. Implications for laboratory collisionless shock experiments are 
discussed. 



Supersonic, well collimated jets are observed from var- 
ious astrophysical objects, such as protostellar disks [![, 
X-ray binaries [2| and active galactic nuclei Q, which 
exhibit wide range of size and radiation power. Despite 
their rich phenomenology, many questions on the physics 
of astrophysical jets, e.g. launching mechanism, collima- 
tion mechanism, role of magnetic field and interaction 
with ambient medium, still remain unanswered. While 
the traditional approaches to address these questions are 
mainly direct observation and theoretical modelling, lab- 
oratory produced jets with proper scaling relations 0, Q 
may provide an alternative platform to study jets on as- 
trophysical scales. In particular, jets formed by laser pro- 
duced plasma have been successfully created by either ir- 
radiating a planar target with one ^concaved" (intensity 




or shining multiple 
101 1 . In both designs 



lower at the center) laser beam [g 
beams onto a cone-shaped target 
ablated plasma plumes are produced at different loca- 
tions on the target and collide on the central vertical axis. 
The radiative collapse effect, which reduces the pressure 
along the axis, then channels the plasma flow along the 
axis and help maintain a jet- like structure [ll| . Here with 
application to laboratory collisionless shock experiments 
in mind, we propose a new design to generate versatile 
laboratory plasma jets with a large dynamic range. It 
utilizes multiple intense high-energy laser beams focused 
onto a simple planar target to form a hollow ring pattern. 

The setup of our numerical simulations is illustrated 
in Fig. [TJ which is based on the Omega laser parameters 
used in recent collisionless shock experiments [l2| . Here 
ten Omega laser beams (of 0.5 kJ each) are used to irra- 
diate a planar plastic (CH) target. Each beam produces 
a supergaussian focal spot at 351 nm wavelength with a 
diameter of 250 /im. The target foil measures 2 mm in 




Side view Top view 

FIG. 1. (Color online) Design of laser-produced plasma jet 
experiment, (a) side view; (b) top view of target surface in 
actual experimental setup (ellipses approximated by circles) ; 
(c) numerical representation in 2D simulation. 



diameter and 0.5 mm in thickness. These beams hit the 
target with an incident angle of 30° from target normal, 
delivering ~ 5 kJ total energy in a 1 ns square pulse. If 
all beams hit the same spot at the center of the target, 
this gives a combined intensity of ~ 10 16 W/cm 2 . We 
carry out axisymmetric radiative hydrodynamic simula- 
tions using the FLASH code [l~3| in a 2D R-Z cylindrical 
domain. This means that even though the actual target 
surface looks like Fig. [TJb) , our idealized numerical sim- 
ulation is really modelling a uniform ring of laser energy 
deposition on the target as in Fig. [Tic) . The radius of the 
ring d is the distance between the center of an individual 
laser beam and the center of the target. The simula- 
tion domain is covered by a uniform grid with resolution 
of (N r x N z ) = (512 x 2560) (i.e. - 2 ^m resolution) 
and open boundaries. The domain is initially filled with 
low density (2 x 10~ 7 g/cm 3 ) Helium to mimic a vacuum 
environment. 

Fig. [2] compares the state of ablated plasma at t = 
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FIG. 2. (Color online) Snapshots of electron density, flow velocity, electron and ion temperatures at t = 1.5 ns for two different 
runs. In one case, all the beams hit the target center (upper panels). In the other case, the focal spot is 800 /im away from the 
target center (bottom panels). 



1.5 ns (0.5 ns after the laser is turned off). In the up- 
per panels, all the laser beams focus on the center of the 
target. In the upper panels, all the laser beams focus 
on the center of the target. As expected, the laser en- 
ergy heats up the target and produces a quasi-spherical 
plasma plume which expands into the ambient low den- 
sity helium. The shock "tip" in the center and two "roll- 
up" s on the side seen in all four sub-panels are proba- 
bly caused by the angle formed between inclined beams 
and the presence of ambient Helium. Note that in ac- 
tual experiments, which take place in a vacuum, there is 
no shear between CH flow and ambient medium, these 
features would be much less pronounced or even absent. 
Other than these features, distributions of plasma prop- 
erties in this case clearly show quasi-spherical uniformity. 
In contrast, beams in the bottom panels form a ring of 
800 /im radius. And they produce a cylindrical plume 
along a ring around the center. As it expands, the plume 
collides on axis and create a highly collimated outflow 
with high density, high velocity and high temperature 



along the axis (see also Fig. [3]). In panels (e), (f) and (h), 
we can see the formed jet structure is on average about 
one beam size in width, and its boundary is marked by 
a sharp transition in plasma conditions. This feature is 
less prominent in electron temperature plot. Compared 
with the case with d = 0, although the laser irradiation 
intensity is lower (the same total energy gets more spread 
out), the jet head travels a similar distance [z ~ 0.4cto) 
because of higher velocity along the axis. The formation 
of the high-density, high-velocity jet is due to the rocket 
nozzle effect and the long, thin, channel visible on the 
z-axis in panel (e) is caused by the rarefication from the 
on-axis plasma collision. 

Evolution of the plasma jet is depicted in Fig. [3]where 
the laser beams from a ring of 800 fim radius. By t = 4 
ns (3 ns past the end of laser energy input), the jet has 
evolved into a structure as wide as the target. And the jet 
head has gone way past the boundary of our simulation 
box. 

To study quantitatively the effects of hollow beam ra- 
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FIG. 3. (Color online) Evolution of the plasma jet produced 
with the radius of on-target beam ring pattern being (upper) 
and 800 /im (lower). 



dius on the jet properties, we plot in Fig.[4}a)-(f) the elec- 
tron, ion number densities, electron, ion temperatures, 
flow velocity and Mach number as a function of time for 
four different runs. The Mach number is computed as 
the ratio of flow velocity to sound speed (c s = y/jP/p). 
All these quantities are taken from 4 mm above the tar- 
get surface and averaged over one beam size (250 pm) 
around the axis. Fig.|U[a) and Fig.HJb) show that as the 
laser beams become more separated, both electron and 
ion number density on the axis increase. All four curves 
top around t = 5 ns but from d = to d = 800 pm 
the maximum density goes up by almost one order of 
magnitude. This is due to the fact that jet collimation 
becomes more prominent with larger beam separation 
(plasma more concentrated on the axis). In Fig. HJc) 
and Fig.^d), particle temperatures are also higher most 
of the time, although the increase is not as drastic as the 
density. The peak velocity of the d = 800 pm case is 
approximately 60% higher than the d — case. How- 
ever the corresponding Mach number is lower because of 
higher ion temperature and higher sound speed. Fig. 0] 
shows that we can produce a much larger dynamic range 
of density, temperature and velocity in a laser-driven out- 
flow, by varying the radius of the hollow laser system. 
This allows a much more versatile and flexible platform 
for all sorts of potential laboratory astrophysics experi- 
ments, from collisionless shocks, to shear flows, jet prop- 
agation and interaction with ambient media. 

The d = benchmark case has actually been realized 
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FIG. 4. (Color online) Average plasma conditions as a func- 
tion of time measured at 4 mm above the target surface along 
the axis. They are the electron, ion number densities (a, 
b), electron, ion temperatures (c, d), flow velocity and Mach 
number (e, f). Different line types represents different degrees 
of beam separations. The bottom two panels show collision 
mean free path (black solid), characteristic electrostatic (blue 
dashed), electromagnetic (red dot-dashed) instability length 
scales and ion skin depth (green dotted) for two different laser 
beams separation cases. 



in recent Omega experiments to study collisionless shocks 
12j, Il4j . There two plasma plumes were produced by ir- 
radiating ten Omega beams each on two facing planar 
targets. The hope was to achieve a laboratory collision- 
less shock from the plume collisions which can be used 
to study astrophysical collisionless shocks. To evaluate 
the viability of shock formation, three interesting length 
scales need to be considered. They are the ion collision 
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and characteristic electromagnetic instability length scale 
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where A z , Z and n z are the ion mass, average charge and 
ion density, v is flow velocity, W(eV) = 5.2 x 10~ 13 A z v 2 
is the kinetic energy per ion, c/uj p i is the ion skin depth 
and T e is the electron temperature. K\ and K-i are two 
poorly constrained numerical factors and we follow Ref. 
[121 ] to take K\ = 30, K 2 = 100. A necessary but in- 
sufficient shock condition is that \ m fp has to be greater 
than both l* ES and l* EM - These length scales, plus the 
ion skin depth at 4mm from the target for cases with 
d = and d = 800 /urn are shown in the bottom panels 
of Fig. |U In the zero radius case, the above condition is 
satisfied from t ~ 2 ns to t ~ 8 ns whereas in the 800 
fj.ui radius case, the condition is satisfied up to t ^6.5 ns. 
Since both electron and ion density peaks around t = 5 
ns, this slight increase in collisionality has a minor ef- 
fect since, we can still take advantage of the density and 
temperature enhancements brought about by the larger 
laser beams ring radius. On the other hand, there are 
two important advantages of the 800 /jm radius case: (a) 
Both the electrostatic and electromagnetic length scales 
are now much smaller in absolute values, implying poten- 
tially thinner shock structures than the 4 mm collision 
distance [14[ . (b) The ion skin depth (bottom green dot- 
ted curve) is also much smaller, which gives the required 
plasma instabilities more room to operate (i.e. easier to 
form collisionless shock). In summary, the higher densi- 
ties and temperatures of colliding jets driven by the 800 
fj,m radius beams are much more favorable for the forma- 
tion and study of collisionless shocks. Even though the 
Mach number is reduced from ~ 18 to ~ 7 (Fig. [H[f)), it 
is still high enough for strong shocks to form. 

Regarding to the "optimal" hollow beam radius for 
jet production, there really is no definitive answer as it 
largely depends on what plasma platform one wants to 
create, and what astrophysical processes one intends to 
study, specifically. Some applications may require higher 
density, whereas high Mach number might be more im- 
portant for others. Plus there are always constraints im- 
posed by the experiment facility. The general trend we 
observed from our 2D simulations is the following: as the 
hollow laser beam radius is increased, the ablated plasma 
plumes collide on axis at a later time, leading to a bet- 
ter collimated, more axially condensed plasma jet, thus 



higher density, higher temperature and higher velocity. 
This is however, countered by the effect that as the laser 
beams get more spread, the laser intensity decreases for a 
given laser energy, leading to lower plasma ablation rate. 
Moreover, when the beams get too far away from each 
other, the plasma plumes would not have even collided 
yet when the laser is turned off. This means that on- 
axis collision would carry less and less energy and the 
resulting plasma jet structure would gradually fade away 
if the beam separation keeps increasing. At least for the 
Omega laser parameters, optimal density and tempera- 
ture increases at 4 mm distance from laser target seem 
to be achieved for hollow beam radius of ~ 800 fim. Be- 
cause of the many competing factors, each laser config- 
uration (e.g. NIF beams) must be studied separately 
for its optimal radius. However, our results demonstrate 
clearly that by varying the hollow laser beam radius, we 
can achieve much larger dynamic ranges for laser-driven 
plasma jets. At least for the collisionless shock experi- 
ment, a higher density allows the jet-jet collision to take 
place at a larger distance from the laser target. This 
increases the space and time scales to allow the shock 
to form and propagate, facilitating the observation and 
diagnostic of the experiment. Another beneficial effect 
is the generation of magnetic fields. Higher density and 
temperature means higher plasma pressure, which helps 
to generate stronger magnetic fields. One caveat of our 
study, though, is that the results so far have all been 
based on 2D simulations which assumed pure asymme- 
try in azimuthal direction. This is not true in reality. We 
are current working on full 3D simulations where the re- 
sults will be more applicable to laboratory experiments. 
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